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Abstract 



We propose a scenario in the Next-to-Minimal Supcrsymmetric Standard Model in which 
its both lightest scalar and pseudoscalar Higgs bosons have masses around 125 GeV. The pseu- 
doscalar can produce an enhanced 77 decay rate due to light higgsino-like charginos in its 
effective one-loop coupling to two photons. In that case, it should contribute to the 77 rate 
measured at the LHC, but not to the WW/ZZ decay modes. Thus a discrepancy should be ob- 
■ served between the rates in these two decay channels through precise measurements. However, 

I/"") \ the pseudoscalar will stay hidden behind the SM-like scalar Higgs boson in the dominant gluon 

fusion production mode and in order for it to be observable the bbh production mode has to be 
j considered. We analyze the constrained NMSSM with non-universal Higgs sector parameters, 

and identify regions in its parameter space where the lightest pseudoscalar with mass around 
j 125 GeV and a measurable 77 rate in the bbh mode can be obtained. 

^ . 1 Introduction 

5_i . 

Since its observation at the LHC in July 2012 [TJ [2], the CMS and ATLAS collaborations have 
accumulated more data and updated their results on the Higgs boson. In the early results, a 
considerable enhancement in the 77 and ZZ rates compared to the Standard Model (SM) prediction 
was noted near ~125 GeV at the ATLAS detector. According to the CMS data, the signal strength, 
ct/ctsM) was consistent with the SM prediction in the ZZ channel but an enhancement in the 77 
channel was observed there as well. However, the figures from both experiments have changed 
in the latest results released after the collection of ~20/fb of data [3j Hj. The signal strengths 
measured by the CMS have now fallen down to SM-like values in the 77 and ZZ decay channels 
with the mean value of the boson mass being 125.6 ± 0.64 GeV. The ATLAS collaboration, on 
the other hand, still reports significant excesses, ct/ctsm = 1-65 ± 0.35 in the 77 channel with the 
mass measurement yielding 126.8 ±0.73 GeV, and ct/ctsm = 1.7 ±0.5 in the ZZ channel with mass 
at 124.3lg4 5 GeV. Moreover, broad peaks consistent with a 125 GeV boson have now also been 
observed in the H — > WW — » 2l2v channel at the two detectors. Importantly, the best-fit signal 
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strength in this channel is very SM-like according to both CMS and ATLAS, with a measured value 
of 0.76 ± 0.21 at the former and of 1.01 ± 0.31 at the latter. 

Since the first announcements of the discovery of the boson, there have been many attempts to 
interpret the observed data in light of various super symmetric (SUSY) extensions of the SM [SHUE! 
El 13 EB HH H3 H3 H31 US] . In the context of the Minimal Supersymmetric Standard Model (MSSM) 
the observed signal can be interpreted as being due to the lightest Higgs boson of the model, h. In 
the MSSM constrained at the grand unification theory (GUT) scale, referred to as the CMSSM [16], 
h can attain a mass around the measured central value only if the SUSY-breaking scale, Mstjsy, 
is much larger than 1 TeV, while also satisfying other important phenomenological constraints. In 
the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [III HE] ( see > e -S-, [13 ED] for a 
review) it has been shown that either of the two lightest CP-even Higgs bosons, hi and hi, can 
easily be SM-like with mass around 125 GeV [211 HI] - In fact in this model it is possible to have 
hi and hi almost degenerate in mass around 125 GeV [9], so that the observed signal is actually a 
superposition of two individual peaks due to each of these, which cannot be independently resolved. 

In the GUT-constrained version of the NMSSM (CNMSSM) [13 E21 EH [24] , in analogy with the 
CMSSM, it has been found that in order to obtain hi as heavy as 125 GeV, MgusY at or above 1 TeV 
is needed even with relevant phenomenological constraints imposed [H]. Alternatively, a SM-like 
hi with mass ~ 125 GeV is easily achievable |14j . Relaxing slightly the universality conditions by 
disunifying the masses of the scalar Higgs doublets mn u and mn d from the scalar mass parameter 
itiq and the soft Higgs trilinear coupling parameters A\ and A K from the unified soft Yukawa 
coupling Aq, makes it relatively easy to obtain SM-like hi or hi around 125 GeV |25j . Here we refer 
to such a model with non-universal Higgs sector parameters as CNMSSM-NUHM. The scenario 
with mass degenerate hi and hi satisfying also other phenomenological constraints has also been 
pursued with much interest in the CNMSSM-NUHM PITO]. 

Even though the latest results from CMS seem to favor a SM-like Higgs boson, those from 
ATLAS do so only partially and it is still possible for the observed boson to be a non-standard 
one. The inconsistencies between the various measurements and fluctuations in the data leave 
ample room for speculation in this regard. Therefore, in this article, we propose a scenario, not 
investigated hitherto, in which the lightest pseudoscalar of the model, a\, is almost degenerate 
in mass with the lightest ~125GeV scalar, h\. Such ai will thus be partially responsible for an 
enhanced signal rate in the 77 channel but, being a pseudoscalar, will not contribute to the WW 
and ZZ channels: 

Y 

i^ 7 (obs) = = Rj n (hi) + R^(ai) ~ 1 + R^{ai) and Rj^ w /zz {ohs) = Rww/zz( h i) - 1 - 

a SM 

(1) 

In the above equation the signal rate is defined as 



, y , u , _ a(Y -» hj) BRjhj ->■ X) 
<r{Y -> h SM ) X BR(h SM X) 



R^ihi) = i^r^ x D „;P (2) 



where X = 77, ZZ/WW and hsM is a SM Higgs boson with the same mass as hi. Y implies the 
various possible Higgs production modes at the LHC, which include gluon fusion (ggh), vector boson 
fusion (VBF), Higgs-strahlung off a vector boson {Wh/Zhj^ and associated production off a heavy 
quark (tth/bbh). Our proposed solution is made possible by the contribution of a light higgsino-like 
chargino, \i > to the effective couplings of ai to a 77 pair, hi is still SM-like in this scenario due 
to a significant singlet component even though tan /3 can take fairly large values [26] . Since this 
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scenario is compatible with a SM-like scalar Higgs boson, it is not in conflict with the recent CMS 
measurements in the ZZ mode \27\ [28] which disfavor the pure pseudoscalar hypothesis. 

We stress that in our scenario the enhancement in the 77 rate due to the a\ contribution will 
not be visible in the ggh production mode. The reason is that the small enhancement achievable in 
BR(hsw-+x) ra ti° wm be negligible in this production mode due to the very small effective coupling 
of the pseudoscalar to gluons which is dominated by the top quark loop. Such an enhancement can, 
therefore, only be noticeable in the bbh production mode which will also get enhanced proportionally 
to the bb decay width of a\. For this reason we shall investigate this particular mode of Higgs 
production here. We emphasize the fact that a measurement of the signal rate due to the bbh 
production channel, which is extremely subdominant for a SM Higgs and is therefore generally 
ignored, becomes enhanced by tan 2 /3 in SUSY models and could lead to a very clear signature of our 
proposed scenario. Furthermore, a difference of the mass measurements in the 77 and ZZ modes 
would also provide a hint for mass degenerate hi and a\. Such a degeneracy would imply that the 
signal observed in the 77 channel should in fact be interpreted as the 'sum' of two individual peaks 
due to h\ and ai, while the peaks in the ZZ/WW modes are because of h\ alone. 

We explore regions of the CNMSSM-NUHM parameter space favoring such a scenario, expecting 
that a discrepancy between 77 and WW/ZZ rates will be seen by CMS and ATLAS collaborations 
with a focussed analysis of the bbh production mode. We investigate the impact of other important 
experimental constraints on the regions of the model parameter space accommodating this scenario. 
They include the limits from direct SUSY searches released by ATLAS with ~20/fb of data as well 
as from the dark matter (DM) relic density measurements. We also require the corresponding 
parameter space to satisfy the recently announced positive BR(B S — > fi + ^~) measurement from 
the LHCb collaboration. 

The article is organized as follows. In Sec. [2] we discuss the possibility of observing an enhance- 
ment in the 77 rate at the LHC due to a ~125 GeV pseudoscalar Higgs. In Sec. [3] we outline the 
model parameter space. In Sec. [4] we describe the experimental constraints applied in our scans, 
present our numerical results and discuss their salient features. We summarize our findings in 
Sec.E 



2 Enhancement in the observed 77 rate due to a light pseudoscalar 

In this section we present some analytical details of our proposed NMSSM scenario in which the 
correlation between the 77 and WW/ZZ rates can be altered. One way to achieve this is with 
mass degenerate lightest doublet-like scalar Higgs, h\, and lightest singlet-like pseudoscalar, a\. 

2.1 The pseudoscalar mass 

We first discuss the conditions that are necessary to obtain a ~125 singlet-like a\ which couples to 
two photons through loops of fermions and charginos only. Starting from the 2x2 pseudoscalar mass 
matrix (after rotating away the Goldstone mode) |19j . one can obtain the approximate expression, 



m 



-3^ SY -^. (3) 



01 ^ M 2 



In the above equation Mp 12 — A(A^ USY — 2ks)v is the off-diagonal entry of the pseudoscalar mass 



matrix, where v = y v u~\~ v d — 174 GeV, with v u and being the vacuum expectation values 
(vevs) of the u-type and d-type Higgs doublets, respectively, and A^J^ Y denoting A x / K at Mstjsy- 
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Mp n ~ [i e gB e ff tan /3, with // e ff = As (s being the vev of the singlet field S), B e g = A^ SY + ks and 
tan/3 = v u /vd, is the diagonal term corresponding to the mass-squared of the doublet-like heavy 
pseudoscalar, ai- The leading term in eq. (|3j) implies that, for positive k, which we will assume 
here, the condition of the positivity of rr? a depends predominantly on the relative signs of /x e ff and 
A K . We explain the effects of negative and positive /i e ff in the following. 

For n c s < (and therefore negative s, assuming positive A), the first term in eq. ([3]) is positive 
if A K > at Msusy- The second term then gives a negative correction to m„ depending on the 
sizes of Mp 12 , which is positive definite, and Mp n , which must be positive for a non-tachyonic 
a-ii requiring in turn B e g < 0. Assuming that the correct m ai is achieved by adjusting the free 
parameters in the leading term, the negative contribution from the second term should be kept 
close to zero. This would require Mp u > Mp l2 . For given tan/3 and /z e ff, 

n is driven by the 

magnitude of B c g, in order to enhance which A\ (which is bounded from above by k\s\) should 
take smaller values. However, A\ at Msusy runs upwards from its GUT value with falling negative 
Aq owing to the contribution from relevant term in its renormalization group equation (RGE) [19]. 
Hence increasing negative Aq diminishes the difference between the two terms in B c g, reducing its 
size and in turn driving Mp n closer to zero. At the same time Mp^, which is a sum of 2k\s\ and 

A\, grows as A\ increases, as opposed to Mp n . Consequently, the ratio . J/ 12 in eq. ([3]) grows with 

' M p,u 

decreasing Aq and, for large negative values of the latter, can result in negative rr? a . Note also 
that the running of A K in turn depends dominantly on ^4^. A K runs upwards with A\ as long as 
the latter is negative. When A\ turns positive A K runs in the opposite direction, owing to its RGE. 
Thus A K in the leading term in eq. ([3j) will have somewhat constrained GUT scale values that can 
yield correct m ai . On the other hand, for fj, e g > 0, the two terms in B e g are both positive and the 
cancellation described above does not occur. 

In summary, the net effect of the interplay between various Higgs sector parameters is that for 
negative n e g the values of Aq at the GUT scale are bounded from below by the condition of the 
physicality of a±. This constraint on Aq causes a slight tension between and m ai , since it is 
well known that to obtain hi which is SM-like with mass ~125 GeV large negative values of Aq are 
required for Mstjsy ~ 1 TeV. For positive n e g there is no such tension because Aq is relatively free 
to take values that give large negative At at Msusy, as long as the correct a\ mass can be achieved 
by adjusting other free parameters. 



2.2 77 decay of the pseudoscalar 

Besides a singlet-like a\ with mass similar to that of the experimentally observed boson, this scenario 
also requires a low mass, m ±, of the lightest chargino. The effective coupling of a pseudoscalar ai, 

Al 

with i = 1,2, to two photons (see, e.g., [221 EQ] ) , is dominated by such a light chargino in the loops 
and can, therefore, be approximated by 

Q C f(77)^ A<?(n), (4) 

X-l 

TV? 

where n = 4 . For r, < 1, which is applicable here, with m ai ~ 126 GeV and the light chargino 

x t 

obeying the lower limit, m x ± > 94 GeV [31 j . the form- factor A"*, 2 (rj) = ^ arcsin 2 yfri |32] in the 
above equation lies in the range 

1 < A\) 2 {n) < 1.2. (5) 
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The coupling of a, to charginos in eq. can be written, following the notation of [19], as 



-^= Pj3 sin 0;7 sin^y — ^j=(Pi2 cos 9u sin fly + Pa sin0[/ cos#y) 
v2 v2 



(6) 



where Ojj, By are the mixing angles for rotating the chargino interaction states to mass eigen- 
states, and Pij are the entries of the mixing matrix that diagonalizes the pseudoscalar mass matrix. 
When the pseudoscalar weak eigenstates Af eak are expressed in the basis {H^i, H u i, Sj) [19] . Pn 
corresponds to H^i, P%2 to H u j and P^ to Si. 

The first term in eq. ([6]) implies that sin^y ~ 1 (yielding a higgsino-like X\)i -Pl3 — 1 an( A 
that larger values of A are needed in order to enhance (77) for the singlet-like a\. Under these 
conditions, the coupling of a\ to bb is naturally suppressed due to the tiny H^j component. On 
the other hand, for the doublet-like pseudoscalar, a2, an enhancement in (77) requires either 
cos 8jj sin By or sin Bjj cos By to be non-negligible. This can be realized only in a very limited region 
of the parameter space where M2 — fj, e g and not too large in order to keep m ± low. Moreover, in 
such a case, the mixing angles in the chargino sector read 



Ouv — arctan 



y/2M w Meff(l + tan/3) 



(7) 



where mw is the mass of boson. The sign of the first term implies that the enhancement can 
only be seen when ai has a leading H^i component so that the term in eq. ([7|) proportional to 
sin 6jj cos 9y is dominant. Evidently, in such a case the a^bb coupling and in turn BR(a2 — > bb) will 
also get enhanced. Consequently, a contribution from 02 will provide no significant excess in the 
77 signal rate, defined in eq. ([2]). 

The above explanation also precludes such a scenario in the MSSM, where the pseudoscalar, 
A, is doublet-like. Besides, as noted in [3H [34] , in the MSSM in order to obtain a lightest CP- 
even Higgs boson, h, with mass around 125 GeV, itia is required to be > 300 GeV, which is 
the so-called decoupling regime of the model. On the other hand, while it is also possible to 
have a ~125GeV H, the heavier CP-even Higgs boson of the MSSM, this can only be achieved 
for 95 GeV < tiia < 110 GeV, in a tiny portion of the 'non-decoupling regime'. This region is, 
moreover, disfavored by the constraints from flavor physics [35^ 136]. 

In the fully constrained version of the NMSSM, unification of A K and Aq at the GUT scale 
introduces tension between the masses of h\ and ai, not allowing both to acquire values <125 GeV 
simultaneously. There, in order to obtain the correct h\ mass, large negative values of Aq are 
necessary so that the mixing term ( m^vsy ~ Msu SY ) can ^ e maxrimze d- A light a±, on the other 
hand requires small A K at MgusY; which in turn implies small A K at the GUT scale, owing to the 
effects of running. Moreover, small values of /i e ff , necessary to obtain light higgsino-like charginos, 
additionally limit the running of A t in the CNMSSM [H]. Therefore, to obtain a SM-like ~ 125 GeV 
h\ and a pseudoscalar with a similar mass and a non-negligible 77 rate one has to look beyond the 
MSSM and the CNMSSM, hence we analyse the CNMSSM-NUHM here. 

Through the mechanism explained above, a more precise measurement of the reduced effective 

C eff (77) 

coupling, C ai (77) = y ? — r> can yield an effective limit on the mass of the lighter chargino 
throughzl 

„ , N , 130 GeV 

C ai ( 7 7)^Ax , (8) 



2 Assuming a singlet-like ai, which implies Pi 3 ~ 1, and a higgsino-like xt so that sin6u,v 
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for m ai ~ 125 GeV. The bound obtained on the mass of \\ is also an effective upper limit on the 
mass of the lightest neutralino, x (= X?)- 

Having described the mechanism for enhancing the 77 decay rate of a±, we now discuss the 
actual quantity used for comparison with the experimentally observed 77 rate. In terms of the 
reduced effective couplings, C 01 (77) and C ai (dd), of a\ to 77 and bb, respectively, the signal rate, 
given in eq. ([2]), can be rewritten, specifically for the aforementioned associated production mode 
with a bb pair in the final state, as 

(aO = ^(^(77) §| - |^| 2 AW/3 (^^) 2 ( f^rg )' W 

where \P[^\ — I ^liv — 2 * s ) v 1 anc j ptotal d ona ting the theoretical value of the total width of a 

SM Higgs boson with the same mass as a±. Note that while larger values of tan/3 are preferred by 
C% 1 in the first equality in the above equation (and also C? (dd) in the corresponding equation 
for i?^(/ii)), there is a further dependence on tan /3 in /T*° which is not as straightforward. 

In the following section we will use eqs. (jH) and ([9]) to obtain an effective upper limit on m ± 
and the mass of x> w-x' m our m °del under consideration. 



3 The CNMSSM-NUHM 

In the fully constrained NMSSM universality conditions are imposed on the dimensionful parameters 
at the GUT scale. This leads to a unified gaugino mass parameter, mi/2, in the model Lagrangian, 
a unified scalar soft SUSY-breaking mass parameter, besides thq and ^4q , to which A\ and A K are 
also unified. These three parameters are thus the only free parameters in the CNMSSM besides the 
dimensionless Higgs coupling A which is taken as an input parameter at Mstjsy, given the correct 
value of the mass of Z boson, rriz- 

In the partially unconstrained version of the model, the CNMSSM-NUHM, the soft masses 
of the Higgs fields, rriH u , m H d and ms, as well as the soft trilinear coupling parameters A\ and 
A K are taken as free parameters at the GUT scale, instead of assuming their unification with rriQ 
and Aq, respectively. Through the minimization conditions of the Higgs potential the three mass 
parameters mn u , mn d and ms at the electroweak scale can be traded for the the dimensionless 
singlet triple-coupling k and the parameters n e g and tan/3. The model is thus defined in terms of 
the following eight continuous input parameters: 

m , mi/2, Aq, tan/3, A, k, /x eff , A x = A K . 

The unification of A\ and A R at the GUT scale assumed here is in general not necessary in the 
CNMSSM-NUHM. In fact, one can argue that the restriction on Aq for [i c r < and the resultant 
tension between and m ai discussed in the previous section can be relaxed by not imposing 
such a condition. In that case, the affect of large A\ can be counter-balanced by increasing A K 
independently, thus still yielding physical a\ solutions. However, this unification condition has 
minimal impact on the allowed parameter space of the model for our purpose, since, as we shall 
see later, we can still exploit the interesting phenomenology of the model while keeping the free 
parameters to a minimum. This is also consistent with the fully constrained version of the model we 
studied earlier [H], where A K and A\ are set equal to Aq at the GUT scale even though ms / tuq. 

In the CNMSSM-NUHM, the upper limit on m, ± and m y , as discussed in the previous section, 

can be obtained from Fig.[U where C ai (77) and i£*,(ai) are shown as functions of m^± in (a) and 
(b), respectively. For all points in the plots we assume 122 GeV < m/i 1)(I1 < 130 GeV. 
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m % ± [GeV] [GeV] 
(a) (b) 



Figure 1: (a) Distribution of points obtained in our scan of the CNMSSM-NUHM parameter space 
in the (m ±,C ai (77)) plane. The dashed line shows the effective upper limit observed in the 
scan. The solid line is based on a perturbative upper limit on A and is shown for comparison, (b) 
Distribution of points in the (m ±, R^L(ai)) plane. The dashed line shows the effective upper limit 
observed in our scan. 



4 Methodology and results 

We perform scans of the parameter space of CNMSSM-NUHM requiring the masses of hi and a\ to 
be 125.5 GeV, with the theoretical and experimental errors taken in quadrature as 3 GeV around 
this mean. We impose the latest 95% confidence level (CL) exclusion limit on SUSY provided by the 
ATLAS collaboration [37] 1^1 We additionally include Gaussian likelihoods for the most significant 
6-physics observables, with their measured mean values and errors taken as: 

• BR (B s m + aO = (3.2±i;| ± 0.32) x 10" 9 , 

• BR (B u tv) = (1.66 ± 0.66 ± 0.38) x 10~ 4 , 

• BR (B -> X s7 ) = (3.43 ± 0.22 ± 0.21) x 10" 4 and 

• AM s s = (17.72 ± 0.04 ± 2.4) ps^ 1 . 

For testing the compatibility of the regions of interest against the direct detection cross section, 
<7p, we use the XENON100 90% CL exclusion limits [38]. Note that we neglect the constraint 
here, since it is well known that the regions where correct can be obtained in the parameter 
spaces of SUSY models with unification of squark and slepton soft masses are strongly disfavored 
by the direct SUSY searches at the LHC [391 HE]- 

The numerical analysis was performed using the BayesFITS package which engages several ex- 
ternal, publicly available tools: MultiNest [10] for sampling of the CNMSSM-NUHM parameter 

3 These limits were originally published assuming the CMSSM. However, it has been verified in [141 17] that they 
have negligible dependence on the Higgs sector parameters and are, therefore, applicable to any R-parity conserving 
SUSY model with unified mo and ra\i%. 
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space; NMSSMTools v3.2.4 [3T] for computing SUSY mass spectrum, Higgs BRs and reduced cou- 
plings, as well as AMb s for a given NMSSM point; Superlso v3.3 02] for calculating BR (B -> X s7 ), 
BR(B S — > and BR(B U — > tv). DM observables such as the relic density and o^ 1 are calcu- 

lated with MicrOMEGAs v2.4.5 03]. 

As noted in the previous section, the scenario under consideration requires low values of 
giving a light higgsino-like xt anc ^ correspondingly a x with significant higgsino component. In fact, 
the CNMSSM-NUHM parameter space conforming to this scenario can be divided into three main 
regions depending on the composition of x : i) the Focus Point (FP) region, ii) the higgsino region 
and iii) the singlino-higgsino region. Below we discuss the results for each of these regions separately. 
We note here that in all our results we show points with neutralino relic density, Q x h 2 , lying in the 
±2cr range, 0.087 < n x h 2 < 0.137, after taking into account 10% error on the theoretical calculation. 
Also, we use slightly extended range of the allowed Higgs mass, 122 GeV < mh 1>ai < 130 GeV, 
compared to the mass measurements of the observed boson at the LHC in order to allow larger 
values of theoretical errors. Finally, for all the points considered, h\ is always SM-like, with 
R^) ~ 1 and R^ z/ww {h x ) ~ 1. 

4.1 The focus-point region 

A light neutralino with mixed bino-higgsino composition can generate correct DM relic density, 
Q x h 2 , in the so called focus-point (FP) region of minimal SUSY models in general [33]. We 
observed in our preliminary scans that this region better satisfies the constraints from XENON100 
and BR(6 — > sj) measurement with /i e g- < 0, so we shall pursue this case here. In Fig.[2^i we show 
the region in the (mo,™^) plane generating a light a\ (122 GeV < m ai < 130 GeV) and x with 
a dominant bino and a small higgsino component. In the figure light blue squares correspond to 
points with 1 < R™ (hi + a x ) < 1.15 and green squares to points with 1.15 < R£L(h\ + a\) < 1.3. 
We see that while large values of mo are favored in order to enhance the mass of h\, is typically 
low, which is necessary for producing a mixed bino-higgsino X- This region, however, lies very close 
to the current 95% CL exclusion limit from ATLAS (also shown in the plot) and should potentially 
be tested soon. 

In Fig. [2b we show the favored ranges of the Aq and tan /3 parameters, tan /3 is almost always 
> 5 to allow enhancement in the h\bb coupling in our considered Higgs production mode, as noted 
earlier. However, we see in the figure that for high positive Aq tan/3 is limited to small values, 
< 15. The reason is that large tan /3 results in an enhanced Yukawa coupling of h\ to bb and rf. 
Consequently A\ runs downwards rapidly from its GUT scale value (we shall see below that large 
positive Aq coincides with negative A\) to more negative values at MgusY- This in turn causes A K 
to run upwards from its GUT scale value, raising m ai beyond its desired value. The effective upper 
bound on tan/3 is relaxed for lower |^4o|> when the running is slower. 

In Fig. [2k we show the distribution of the parameters Aq and A K = A\ at the GUT scale. We 
notice that Aq stops at much smaller negative values than it would be expected to take in order 
to maximize . This is due to the same reason as explained in Sec. [2j The main contribution to 
m ai comes from the leading term in eq. ([3|). Since in this region we assume /j. e s < 0, Aq is strongly 
bounded from below in order to minimize the affect of the second term there. Such negative Aq, 
by causing positive A\ to run upwards, also pushes A K to somewhat large positive values at the 
GUT scale resulting in its small positive values at MsusYj since it runs in the opposite direction to 
A\. This is also the reason why no points are visible in the region with negative A K and negative 
Aq, but one can see some points with positive A K and positive Aq. Finally, for negative A K large 
positive Aq can be reached, since such values of Aq drive positive A\ at the GUT scale downwards, 
which in turn causes A K to run upwards to positive values at Mstjsy- 
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Fig.[2]i shows the ranges of A and k corresponding this scenario. While k is widely distributed, 
the allowed range of A is subject to a three-way tension. Large A is favored in order to obtain an 
enhancement in the coupling of a\ to X\ but the condition to obtain a SM-like hi, on the other 
hand, prefers smaller values. The small-to-intermediate range of A seen in the figure is then a result 
of the compromise between these two conditions and, additionally, of the requirement to achieve 
the desired m ai by generating s (= // e ff/A) of the correct size. 

In Fig. [2^ we show how the points in the FP region yielding such a scenario fare against 
the XENON100 limits. Also shown in the figure are the actual 90% CL exclusion limits from 
XENON100 as well as the 90% CL limits expected from the LUX 05] and XENON1T 06] exper- 
iments. Note that m x is bounded from below by the ATLAS limit on miji in this region as it is 
bino-dominated. We see that a majority of the allowed points with an enhanced 77 rate lie below 
the XENON100 line. Most of this region, however, lies above the LUX limit, while the XENON1T 
data should be able to test almost all of it. Fig.[2f shows the allowed masses of h\ and a\ in this 
region. We see that m/ ll is always lighter than 124 GeV, which is a consequence of the not very 
large values of negative Aq allowed in this region, as discussed above, ai, on the other hand, can 
easily have a mass around 125 GeV. 

Overall, we notice only slight enhancement, up to ~25% in the 77 rate compared to the SM 
expectation in this region of the CNMSSM-NUHM parameter space. The reason being that m ± 
is not allowed to take small enough value due to the lower bound on the mass of x discussed above 
(m x ± ~ fj, B g > m x ). BR(B S — > /i + /i~) in this region varies between 2 x 10~ 9 and 5.5 x 10~ 9 , which 

is within 2a of the experimentally measured value 3.2 x 10 -9 , taking into account the theoretical 
error (as in [14j). On the other hand, BR(fo — > 57) takes values between 3.1 x 10~ 4 and 3.7 x 10~ 4 
and hence is always close to the experimental value. This region, owing mainly to the facts that 
m^j finds it difficult to reach the experimentally observed value and that the reduced 77 rate of oq 
barely exceeds 0.25, is the least favored of the three regions explored here. 

4.2 The higgsino region 

A nearly pure higgsino- like neutralino can generate large enough £l x h 2 only if m x ~ /x e g ~ 1 TeV [7], 
but such high values of // e fr will not yield the desired enhancement in the a± — > 77 rate. Therefore, 
in order to obtain a sizeable enhancement one has to relax the condition on neutralino relic density 
(thereby allowing low fi c g and, therefore, Q x h 2 to be too low). One can assume that a neutralino 
contributes only partially to the relic abundance of the universe beside some other DM candidate. 
In that case Q x h 2 = £^ tota i/i 2 , where £ is the fraction of the total relic abundance produced by x 
and Ototai^ 2 = 0.112. Another possibility is that the entire relic abundance is due to an alternative 
DM candidate particle in the model. Often considered examples of such an additional/alternative 
DM candidate are gravitino (see, e.g., [U] for recent analyses in the MSSM) and/or axino [48j . 
The first (second) of these candidates is (not) tightly constrained by Big Bang Nucleosynthesis but 
both are likely to be allowed in this region due to the low neutralino yield at freeze-out. 

In Fig. [3^, we show the region in the (w<o, m^) plane generating a light ai and a higgsino-like 
X- Light blue squares correspond to points with 1 < R^ (hi + ai) < 1.15, green squares to points 
with 1.15 < R™y(hi + ai) < 1.3, red squares to points with 1.3 < R®L(hi + a\) < 1.45 and green 
squares to points with 1.45 < R^ (hi + ai). Large values of mo are preferred in order to enhance 
m/ ll through radiative corrections from the SUSY sector, myj also takes large values in order to 
minimize the bino component of x with large lJ- e s~' m i/2 splitting. 

In Fig. [3b the favored ranges of tan /? and Aq parameters are shown. We see that the en- 
hancement in R 1 ^ (hi + ai) decreases as tan/3 increases. The reason for this is as follows. The 
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and m ai obtained in this region. See text for details. 
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enhancement in i?™,(ai) grows with A, according to eq. 0. However, large values of A can only 
give correct for not too large values of tan f3. This is because larger values of tan /3 result in an 
enhanced Yukawa coupling of h\ to bb and rf. This will make A\ run upwards faster from its GUT 
scale value, which in turn causes A K to run downwards to larger negative values. That will result 
in a decrease in m^, since it has a significant singlet component, while m ai increases. We point 
out here that an opposite affect of tan f3 was noted in the FP region due to the fact that there Aq 
was positive which made A\ and A K run in the opposite directions to those here. 

Since for this region /i e g > 0, the effective constraint on Aq from A\ applicable in the FP region 
is relaxed here and larger values of negative Aq are allowed so that m/ 11 can be maximized. The 
interdependence of Aq and A\ = A K is further illustrated by Fig. [3b- A K can take much bigger 
positive values at the GUT scale, even though it ought to be negative at Msusy- Moreover, a 
considerable number of points is visible for negative A K and negative Aq up to ~— 2 TeV in contrast 
with the FP region. The reason is, again, the fact that here negative Aq makes A\ run upwards 
between the GUT scale and Msusy and since the latter is negative it also drives A K upwards to 
smaller negative values. Naturally then, negative Aq should not be too large or A K at Msusy will 
be driven positive. For positive A K only large negative Aq values are allowed because A\ is also 
positive now and so the running of A K switches to downward direction. Positive Aq solutions are 
also possible, as in the FP region, as long as they don't yield positive A K at Msusy- 

In Fig.[3]i the ranges of the paramaters A and k favored by our proposed scenario are shown 
for this region. While k is as widely distributed as in the FP region, comparatively slightly larger 
values of A are possible in this region since s is more free to vary owing to the less constrained A K . 
In Fig. [3k we show the distribution of the points in the {m x ,^a^) plane for this regionQ Almost 
all the points obtained in this region lie below the XENON100 line and a portion of these points 
even lies below the projected XENON1T line. Since x 1S almost purely higgsino here, the smaller 
its mass the bigger the enhancement in F& (a±) is generated. Moreover, this region corresponds 
to large values of tuq and mj^ so that the squarks and gluinos are always much heavier than 
the current LHC reach. Nevertheless, as discussed earlier, a more precise measurement of ic* 
could still introduce limits on m x and rn x ±. Such derived upper limits are, therefore, especially 
interesting from the experimental point of view. 

Finally, as a result of larger allowed values of negative Aq, h\ as heavy as 129 GeV can be 
obtained in this region, as can be seen in Fig.[3]f. a\ mass is evenly distributed in the defined range, 
almost always showing a large enhancement in the 77 rate. Finally, compared to the FP region, a 
majority of points in this region show a big enhancement, up to ~50% above the SM expectation, 
in the 77 rate in this region. We also note here that both BR (B s — > fi + fi~) and BR(6 —> 57) always 
lie around their respective SM values. 

4.3 The singlino-higgsino region 

This region is defined by x being an admixture of a large higgsino component and a smaller 
but important singlino component. Owing to the significant singlino component (20%-30%) the 
neutralino will interact very weakly with matter and will thus result in large relic abundance 
unless it has a small mass and consequently large annihilation cross-section. In Fig.Ek we see the 
distribution of this region in the (?7io> m i/2) plane. The color assignment of the points is the same 
as in the higgsino region. Large values of m,Q and intermediate values of m 1 / 2 are favored, again, 
for allowing a higgsino- like neutralino with smaller positive fi e s. In Fig.Sb we see that once again 

4 The figure assumes that \ is on ly responsible for a small portion of the observed relic abundance. For the points 
obtained in the scan £ < 0.05. 
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tan/3 spans a fairly wide range but its small to intermediate values are favored by R!*L(&i), for the 
same reasons as in the higgsino region. Aq almost always takes large negative values, in order to 
maximize m/ ll . In contrast with the higgsino region, there are hardly any points corresponding to 
positive Aq. 

In Fig. St w e show the distribution of points in the (Aq,A k ) plane and in Fig. HE in the (A, k) 
plane. The observed range of k is smaller compared to the previous regions while that of A is 
larger. This is in fact the main feature distinguishing this region from the higgsino region in terms 
of the parameter space of the model. Larger values of A and small values of k are required to 
maximize the singlino component of as the 5x5 term in the neutralino mass matrix is equal to 
Iks (= 2K/U e ff/A). Hence the smallness in k has to be compensated by slightly larger values of A K 
compared to the previous two regions, as noted in Fig. [lb, for obtaining the correct m ai . Finally, 
like the other two regions, the largest enhancement in I& [a\) is obtained for largest possible values 
of A. 

In Fig.0fe we show the (m x , cr^ 1 ) plane for this region. A large number of points satisfying the 
XENON100 limit lies below the projected 90% CL XENON1T limit. Note also that since much 
smaller m x and consequently m^± is favored by this region compared to the other two, almost all 

the points below the XENON100 line have a highly enhanced B^{a\), since Xi a ^ so appears in the 
denominator of eq. ([9]). This is also the reason why such points are achievable even with relatively 
small values of tan /3, as seen in Fig.Ub earlier. This region is thus the one yielding maximum 
enhancement, up to ~60% or so, in R™L(hi + a\) out of the three regions discussed here and is, 
therefore, the most favorable fo all. 

In Fig.Hf we show the distribution of versus that of m ai . We note that this region yields 
mostly larger compared to the previous two regions, which is due to the combined effects 
of large negative Aq as well as larger allowed values of A. Another distinguishing feature of this 
region is that hi can also be almost mass degenerate with h\ and a±, implying in that case a 
'triple degeneracy' among the Higgses. Again, while mass degenerate hi and hi can explain the 
enhanced 77 rate in the ggh production mode in the ATLAS data, in order to test the additional 
degeneracy with a\ one will have to explore the associated production mode of Higgs with bb. In 
the bbh production mode, such hi can further contribute ~20% of the measured 77 rate. Finally, 
BR(B S — > in this region varies between 3 x 10~ 9 and 3.8 x 10~ 9 while BR(6 — > s"f) lies in 

the 2.8 x 1CT 4 to 3.3 x 1CT 4 range. 

To conclude, in Fig.[5ti we show the range of m x across all the regions for which an enhancement 
in B^{h\ + ax) was obtained in our CNMSSM-NUHM scan, and its compatibility with the current 
and expected limits on a^ 1 . These regions are identified seperately in Fig. [5b, again, in the (m x , ^Up 1 ) 
plane, where the orange squares denote the FP region, dark green squares the higgsino region and 
brown squares the singlino-higgsino region. 

5 Summary 

We have proposed a scenario in the NMSSM in which the lightest pseudoscalar of the model, a\, 
as well a SM-like scalar Higgs boson, hi, both have masses around ~125GeV. The pseudoscalar 
could be distinguishable from the scalar at the LHC in the associated production mode with a bb 
pair in the final state. This is because it will contribute significantly to the observed signal rate 
in the 77 channel but, since a pseudoscalar doesn't couple to W and Z bosons, the measured rate 
in the WW/ZZ channels will be due only to the scalar and therefore SM-like. The observable 
enhancement in the 77 rate is made possible by a light chargino entering in the one-loop effective 
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coupling of a\ to two photons. We have discussed the conditions necessary to obtain a\ with the 
correct mass as well as an enhanced 77 rate. We have argued that due to very specific requirements 
on the composition of a\, which should be singlet-like, and of the light chargino, which should be 
almost purely higgsino-like, such a scenario cannot be realized in the MSSM and is extremely 
unlikely in the fully constrained NMSSM. 

We have, therefore, analyzed the CNMSSM with the universality conditions lifted in the Higgs 
sector to study our proposed scenario. We have scanned the parameter space of this model in order 
to look for regions that can allow both X\ an d cl\ with the desired masses and compositions. We 
have found that these regions can be divided into three broad cases, based on the composition of 
the neutralino which, owing to the condition on X\, should also have a large higgsino component. 
These regions include the FP region, where x is a bino-higgsino mixture, the higgsino region, where 
it is almost purely higgsino, and the singlino region, where it is higgsino-dominated but with an 
admixture of the singlino. The region least favored by enhancement in the 77 rate of a\ is the FP 
region, where it only reaches up to ~25%, while the most favored one is the singlino region, where 
the enhancement can be as high as ~60%. 

We noted, however, that such oq is likely to remain invisible at the LHC in the gluon- fusion 
production channel. The reason is that while the effective coupling of a± to the 77 pair gets 
enhanced, the effective coupling to two gluons is still highly suppressed compared to h\. We 
emphasize that a more focussed analysis of the associated Higgs production mode with bb pair, 
which is the least favorable production mode for a SM-like Higgs, is essential. By revealing such a 
pseudoscalar, this production mode could provide a clear signature of our proposed model scenario, 
in particular, and of beyond the SM physics, in general. 
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